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^ . ABSTRACT 

^ ' We study the effect of primordial non-Gaussianity on the development of large-scale 

cosmic structure using high-resolution iV-body simulations. In particular, we focus 
on the topological properties of the "cosmic web" , quantitatively characterized by 
(N ■ the Minkowski Functionals, for models with quadratic non-linearities with different 

values of the usual non-Gaussianity parameter /nl- In the weakly non- linear regime 
(the amplitude of mass density fluctuations (Jq < 0.1), we find that analytic formulae 
derived from perturbation theory agree with the numerical results within a few percent 
of the amplitude of each MF when |/nl| < 1000. In the non-linear regime, the detailed 
behavior of the MPs as functions of threshold density deviates more strongly from the 
analytical curves, while the overall amplitude of the primordial non-Gaussian effect 
remains comparable to the perturbative prediction. When smaller-scale information is 
included, the influence of primordial non-Gaussianity becomes increasingly significant 
statistically due to decreasing sample variance. We find that the effect of the primordial 
non-Gaussianity with |/nl| = 50 is comparable to the sample variance of mass density 
k> ■ fields with a volume of 0.125(/i^^Gpc)'^ when they are smoothed by Gaussian filter 

, at a scale of 5ft.~^Mpc. The detectability of this effect in actual galaxy surveys will 

' strongly depend upon residual uncertainties in cosmological parameters and galaxy 

biasing. 

Key words: Cosmology: early Universe - large-scale structure of Universe - methods: 
N-body simulations - analytical - statistical 
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1 INTRODUCTION 

According to the standard scenarios for the formation of 
large-structure in the Universe, the present-day cosmic den- 
sity field evolves from small-amplitude initial fluctuations 
which are described by Gaussian statistics. The hypothesis 
of primordial Gaussianity is supported by present observa- 
tions of the Cosmic Microwave Background (CMB), particu- 
larly those from the Wilkinson Microwave Anisotrop y Probe 
(WMAP) jKomatsu et al.ll2003l : 'Spcrgcl e t al.ll2007^ . These 
results are consistent with an inflationary origin for the pri- 
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mordial perturbations, since the simplest forms of cosmic 
inflation produce nearly Gaussian fluctuations. 

In order to understand the early Universe in more detail, 
however, it is necessary to measure (or at least constrain) 
the departures from non-Gaussianity that inevitably arise 
at some level during the inflationary epoch. For example, 
the simplest slowly-rolling single field inflation model pre- 
dicts very small levels of primordial non-Gaussianity, while 
multi-field infiation models and models with a non-standard 
kinetic term for the infiation may yield larger effects which 
could be d etected in ongoing o r next-generation observa- 
tions (e.g. iBartolo et al.l 20021 : Bernardeau fc UzanI 20021: 
Lvth et al.l l2003l: iDvaU et all 120041: lArkami-Hamed et al. 



2004l : lAlishahiha et al.ll2004l : lBartolo et al.ll2004l : IChen et al 
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\200i : iBattefeld fc Battefeldll2007l ) . Only when such phenom- 
ena are detected will it be possible to distinguish between 
the hundreds of currently viable variations on the theme of 
inflation by understanding the dynamical behavior of the 
inflation field. 

In order to model the primordial non-Gaussianity 
that might arise during inflation, the following simple 
form including quadratic corrections to the curvature per- 
turbation <& (jBardeen if98Q') during the matter era has 
been often adopted jGangui et al]ll994l : IVerde et"^ ] |2000l : 
iKomatsu fc SpergellbOOlf ii 



' + /nl( 



(1) 



where <j!> represents an auxiliary random-Gaussian field and 
/nl characterizes the amplitude of a quadratic correction 
to the curvature perturbations in a dimensionless way. In 
principle, /nl could be scale-dependent, but current obser- 
vations are not sufficiently sensitive to detect any such vari- 
ation, so a constant /nl remains a useful parametrization 
of the level of non-Gaussianity. Recent analyses of the an- 
gular bispectrum for WMAP provides strong constraints on 
/nl to lie in the range from —54 to 1 1 4 at the 95 percent 
confidence level (iKoinatsu et al.l l2003l : ISpergel et al] l2007l : 
ICreminelh et al.ll2006l l. 

The Large-Scale Structure (LSS) of the distribution 
of galaxies in the Universe provides another poten- 
tially powerful probe of primordial non-Gau ss ianity 
jFrv fc ScherrerdI Il994l: IChodorowski fc Bouched Il996l: 



Verde et al' '2000'; Scoccimarro et alj |2004| : iHikage et all 



2006; Scfusatti & Komatsu 200 J? The three-dimensional 
spatial information arising from LSS is potentially a richer 
source information about primordial non-Gaussianity than 
the two-dimensional information arising from the CMB. For 
example, constraints from upcoming cluster surveys should 
be comparable with current CMB limits and those from 
galaxy surveys, which could be as tight as I/nl] ~ 10 for 
the planned surveys and |/nl| ^0.2 for an all-sky survey of 
galaxies up to re dshift z = 5 (jSefusatti fc Komatsul [20071 : 
IPalal et al.ll2007l ). A variety of large-scale projects of LSS 
observation covering Gpc^ volumes are being proposed, 
such as an extension of the Sloan Digital Sky Survey; APO- 
LSS survey; The Hobby-Eber ly Dark Energy Experiment 
(HETDEX) (|HiU et all |2004|): Wide-Field Mu lti-Ob iect 
Spectrograph ( WFMOS) (|Glazebrook et al.ll2005l): and th e 
Cosmic Inflation Probe (CIP) mission ijMelnick et al.ll2004l l. 
It is consequently important to study the optimal way to 
extract information about primordial non-Gaussianity from 
such surveys. 

The statistical analysis of non-Gaussianity has 
been mainly performed through the calculation of the 
bispectrum (|Verde et al.l l2000l : IScoccimarro et al ] l2004l : 
ISefusatti fc Komatsd 20071 ) . Strong motivation for this is 
that the bispectrum is the simplest statistic al function that 
can measure quadratic non-linearity (e.g. IWatts fc ColesI 
l2003h . Although the quadratic model provides an extremely 
useful benchmark for statistical analysis techniques, one 
must always bear in mind that there are many different ways 
for a random field to be non-Gaussian. In general, there is no 
one statistic that completely characterizes the statistical na- 
ture of a non-Gaussian random field, so a battery of higher- 
order statistics must be deployed. In particular, when the 
full nature of non-Gaussianity is virtually unknown, such as 



is really the case for primordial perturbations, the theoreti- 
cal model assumed should be validated before its parameters 
are constrained. Different statistics reflect different aspects 
of non-Gaussianity so the use of different statistics plays a 
vital role in this kind of consistency check. 

In this paper we use a set of invariant characteris- 
tics of the topology of the cosmic web, known as the 
Minkowski Functionals (MPs). These have already been 
used to describe the morphological proper ties of cos- 
mic d ensity fields in a variety of c ontexts (iMecke et al. 
I994I; 'Schmalzing fc B uchertI Il997l : ISchmalzing fc Gorski 
1998; Hikage ct al. 2003i). Four MPs are defined in three- 
dimensional density fields such as LSS: the volume fraction 
(Vo); surface area (Vi); mean curvature (V2); and Euler char- 
acteristic (V3). 

Using a perturbative approach. iHikage et al.l |200^ de- 
rived analytical formulae for the behavior of the MPs for 
LSS including primordial non-Gaussianity (as a function 
of /nl as given in equation [1]), in addition to the non- 
Gaussianity due to non-linear gravity and galaxy biasing. 
The validity of the perturbative analysis is, however, lim- 
ited to the weakly non-linear regime. Smaller-scale modes 
also contain rich information about the primordial density 
fields, and this could help place more stringent constraints 
on primordial non-Gaussianity. In this paper, we use high- 
resolution A''-body simulations to study the effect of pri- 
mordial non-Gaussianity on the MPs from the mildly to 
strongly non-linear regime. There are two reasons for us- 
ing the full numerical analysis: one is to see how well the 
perturbative formulae describe the simulated MPs to check 
their applicability; the other is to study how the primor- 
dial non-Gaussian effect behaves in the strongly non-linear 
regime and thus to estimate the significance of the effect on 
the MFs. 

The paper is organized as follows. In Section O we re- 
view the perturbative formulae for the MPs. The details of 
the Ai'-body simulations and the computing method of the 
MPs are summarized in Section [3] In Section O we compare 
the perturbative formulae of MPs with simulated results 
to study the primordial non-Gaussian effect in non-linear 
regime. Section [S] is devoted to the summary and conclu- 



2 PERTURBATION THEORY 

We define the MPs of density fields for a given threshold 
f = 5/ao, where S is the density fiuctuation, which has zero 
mean, and ao = (3^)^^^ is its standard deviation. The fc-th 
MP Vki^) can be written separately with the amplitude Ak 
and the function of Vk{v), as 



Vk{iy) = AkVk{iy). 



(2) 



The amplitude part Ak, which depends only on the power 
spectrum P{k, z) of the 3-dimensional fiuctuation field 5 at 
redshift z, is given by 



Ak = 



(3) 

where LJk = n''^^ /r{k/2 + 1) gives ujq — 1, uji — 2, ijJ2 = tt, 
and LOS — 4tv/3. The quantity characterizes the variance 
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of fluctuating fields for i = and that of their derivatives 
for i = 1 given by 



(Ti (2) = 



k^dk 
27r2 



k'^'P{k,z)W'^{kR), 



(4) 



where W represents a smoothing kernel. Throughout the 
paper, we adopt a Gaussian kernel = exp[— (fci?)^] where 
R r epresents the smoo thing scale. 

iMatsubaral (|2003l ) derives the second-order perturbative 
formulae of the MFs using the multivariate Edgeworth ex- 
pansion. According to the formulae, the function Vk{v) is 

(G) 

written with the Gaussian part v^, and the leading part of 
the non-Gaussian term Avk as 



e-'-'/^Hk-iiu), 



Avk{iy) = 
+ 



2/2 



ls^°^Hk+2{u) + ^S^'^Hk{u) 
6 3 



t^S^'^Hk..{.) 



O"0, 



(5) 
(6) 

(7) 



where H„{i') denote the Hermite polynomials. The leading- 
order non-Gaussian term Avkif) is calculated when the 
three "skewness parameters" S^''^ are given. 

The three skewness parameters S^'^li = 0, 1 and 2) are 
computed by integrating the bispectrum B(k\, k2, fcs, z) over 
fci, fc2, and n = (ki ■k2)/(fcifc2) with appropriate weights as 
jHikage et al.ll2006l ) 



dko 



, , , , , dki 
B{ki,k2, ki2,z)W{kiR)W(k2R)W(ki2R), 

fc?fc|(fci + k2+ fj,kik2)B{ki,k2, fci2, z) 
xW{kiR)Wik2R)W{ki2R), 
3 



(8) 



(9) 



16n^af{z) 



dki / dk2 
Jo 

feifc|(l — ^/)B{ki,k2,ki2,z) 

xW{kiR)W{k2R)Wiki2R), 



dfi 



(10) 



where fci2 = |ki + k2| = {kj + k'i + 2^fcife2)^/^. 

Throughout this paper, we neglect the non-Gaussianity 
arising from the non-linearity in relationship between galaxy 
counts and mass (i.e. galaxy biasing) so as to keep the anal- 
ysis as simple as possible. The bispectrum B for the matter 
density fluctuation is then given by 

B{ki,k2, fcs, z) = -Bpri(fci, fc2, fca, z) + Bgr{ki,k2, fes, z), (11) 

where Bpri and Bgr represent the contributions from pri- 
mordial non-Gaussianity and non-linearity in gravitational 
clustering respectively: 



Bpri(A:i, fc2, fcs, z) 
Bgr{ki,k2, kg, z) 



P{ki,z)P{k2,z)M{k3) 

M{ki)M{k2) 



- D{z) 

+ (eye.)] , 

= 2[F2(ki,k2)P(fcl,z)P(fe2,2) 

+ (eye.)] , 



(12) 
(13) 



where D{z) is the growth rate of linear density fluctuations 
normalized such that D{z) 1/(1 -f- z) during the matter 
era. The function M{k) and F2(ki, k2) are time- independent 
kernels describing mode-coupling due to non-linear cluster- 
ing of matter density fluctuations in the weakly non-linear 
regime. These are given by 



M(fc)=2fc^^ 



F2(kl,k2 



5 ki k2 

7 2fcifc2 



fki k2\ 
\k2 kJ 



+ 



2 (kl ■ k2) 

T ^\ ^2 



(14) 



(15) 



We adopt the linear transfer function T{k) by 
lEisenstein fc Hul l|T99al . In comparison with numerical 
simulations, we use the power spectrum of the simulations 
(the details are explained in the next section) at z* — 76.97 
for a theoretical input of the power spectrum P{k,z*) and 
then give the power spectrum at z as 



P{k,z) = 



D^{z*) 



P{k,z') 



(16) 



3 METHODOLOGY 

3.1 Numerical Simulations with Primordial 
Non-Gaussianity 

The A'^-body simulations with primordial non-Gaussianity 
that we us e for this analysis are those described in 
iGrossi et al.l |20o3). These simulations employ 800"^ dark 
matter particles in a periodic cubic box with a side length 
of 0.5/i~^Gpc. The cosmology of our simulations is a flat 
ACDM model with mass density parameter ilm = 0.3, 
baryon density parameter fib = 0.04, Hubble parameter 
h = 0.7, primordial power-law index = 1, and erg = 0.9. 

The initial particles are perturbed from an initially ho- 
mogeneous "glass- like" distribution. The primordial non- 
Gaussianity is incorporated into a Gaussian-random field 
with the ab ove cosmology in the form of equation ([TJ. 
iGrossi et al. I (2007 ) explored 7 different scenarios with 
/nl = 0, ±100, ±500 and ± 1000. We have analyzed aU of 
these simulations, but for brevity in this paper we only 
present results for the Gaussian simulation with /nl = and 
the two extreme non-Gaussian cases /nl = ±1000; results 
for the other simulations with /nl = ±100 are intermediate, 
as expected. 

After Fourier-transforming the primordial non- 
Gaussian field, the dark matter particles are displaced on 
the initial grid assuming the Zel'dovich approximation. 
The simulations are started at z « 100 and the subsequent 
gravit ational evolution is simulated with the GADGET-2 
code (|Springel et al.ll2005l ). The Triangular-Shaped Cloud 
method is used to assign densities onto 512^ grids. After 
Fourier-transforming the grid data, we multiply by the 
Gaussian kernel exp[— (fci?)^], and then transform them 
back to real space. 

It is instructive first to examine the visual morphology 
of the clustering pattern. Fig. [1] shows maps of slices of the 
mass density field with /nl = (middle-row panels) and the 
relative residuals between /nl = ±1000 and /nl = (left 
and right panels). The residual for the map with /nl = 2;, 
Apx, is calculated at each pixel as 
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fnl = - 1000 fill = fnl = +1000 




Figure 1. Slice maps of simulated mass density fields a.t z = 5.15 (top), z = 2.13 (middle) and z = (bottom). The number of pixels 
at a side length is 512 (500/i~^Mpc) and that of the thickness is 32 (31.25/i~^Mpc). The panels in the middle row show the log of the 
projected density smoothed with a Gaussian filter of 10 pixels width, corresponding to 9.8/i~^Mpc. The left and right panels are the 
relative residuals for the /NL=ilOOO runs (equation 1171 ). Each panel has the corresponding color bar and the range considered are 
different from panel to panel. 
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= (p:c - Po)/po (17) 

where px is the number density of mass particles for the 
map with /nl = x. The field is smoothed with a Gaussian 
filter 10 pixels wide (i.e. 9.8/i"^Mpc). The redshifts of the 
maps are 5.15, 2.13 and from top to bottom respectively. 
Similar density structures in the mass distribution appear 
in the residual maps with their contrast at same (inverse) 
sign for positive (negative) values of /nl. For example, a 
large void structure at the right-center in the density map 
also appears in the residual maps. This is because the higher 
density region is initially more (less) enhanced in the positive 
(negative) /nl , as predicted by the local model of primordial 
non-Gaussianity in equation ([T]). 



3.2 Computation of Minkowski Functionals 

The computational method we use for calculating MFs of 
data defined on a grid is based on ideas from integral geom- 
etry, rather than the alternative more cumbersome approach 
of using the differential properties of bounding surfaces. In 
our case the calculation reduces to counting the numbers 
of vertices, edges and sides of the elementary cells cover- 
ing t he structure (|Goles et al.lll996l : ISchmalzing fc BuchertI 
1 19971 ). The range of v is from —3.6 to 3.6 with an equal 
binning width of 0.2. The MFs measured from numerical 
simulations often deviate from analytical predictions even 
for Gaussian realizations due to subtle pixeliz ation effects. 
However, as pointed out bv lHikage et al.l (|2006l ). pixelization 
effects become negligible when computing the difference be- 
tween Gaussian and non-Gaussian MFs. Therefore we focus 
on /S.Vkii'i) (i denoting the binning number of u) that we 
compute as follows: 

(i) We compute the MFs for non-Gaussian simulation 
data Vk and then divide them by their amplitudes Ak 
(equation [3]) to obtain normalized MFs v^- The ao and 
(Ti in Ak are computed from the density fields of the 
simulations. 

(ii) The MFs for Gaussian fields are computed in the 
same way and then divided by their amplitudes A^ where 
the values of cro and o\ are computed from each realization. 
The same cosmological parameters as the A'^-body simula- 
tions are adopted. The normalized MFs v), are estimated 
by averaging MFs over 10 Gaussian realizations. 

(iii) The difference ratio Aii^ is computed by 
Avk^Vk-vfK (18) 



4.1 Agreement with perturbative formulae in the 
weakly non-linear regime 

Fig. [2] shows examples of MFs Vk (left panels) and the dif- 
ference ratio Aufc (right panels) for simulated mass distribu- 
tions in the weakly non-linear regime. We smooth on a scale 
R = 10/i~^Mpc which, at z = 3.96, marks the transition 
to the non-linear regime since the variance of the smoothed 
density fluctuation cro — 0.1. The different symbols show the 
different /nl of and ±1000. The error-bars represent the 
sample variance estimated from 1000 Gaussian realizations 
with the same R, z and box-size as the simulations. The 
perturbative formulae discussed above are plotted with lines 
for comparison. Results for the simulations with /nl = ±100 
and ±500 are found to be linearly scaled between those with 
/nl = and ±1000. 

The theoretical curves reproduce the features of the 
simulated MFs very well. We quantitatively estimate the 
agreement between the simulation results l\v''^™\ui) and 
the perturbative formulae Au^^^' (i^i) by calculating the 
root-mean-square (rms) differences averaged over i. Table [T] 
lists the differences for each MF at different redshifts z (but 
R is fixed to be 10/i~^Mpc). The differences are less than 
a few percent relative to the amplitude of each MF (equa- 
tion [3]) when cro < 0.1 and remains at the 10-percent level 
when ctq ~ 0.2. We also estimate the rms differences divided 
by the rms of At;^^'^'(i/i) averaged over i. These quantities 
represent the extent to which the theoretical predictions im- 
prove going from linear theory to (2nd-order) perturbation 
theory. The differences between the 2nd-order perturbative 
predictions and the numerical simulations is 0.15 ~ 0.41 
times smaller than those corresponding to linear theory at 
cro < 0.1. These results are cons istent with the previous 
analysis bv iNakagami et al] (|2004 ). 

The differences between theory and simulations are 
quite small compared to the sample variance. However, there 
is a systematic feature, seen in the asymmetry of Vo and V2 
with respect to = 0; the perturbative predictions are sym- 
metric. There are three possible explanations for this effect. 
One is that higher-order contributions - i.e. beyond 2nd- 
order - are significant. Another possibility arises from the 
use of the Zel'dovich approximation to set the initial con- 
ditions of the simulations, which may be responsible for an 
extra contribution to higher-orde r statistical propert ies of 
clustering arising from transients l|Crocce et al.l [20061 ). The 
other reason is the fact that the multivariate Edgeworth 
expansion which is the basis of perturbation formulae has a 
limite d range of validity, especially a t values of v larger than 
unity (iBernardeau fc Kofmanlll995l '). These effects must be 
considered carefully when comparing with real survey re- 
sults. 



4 RESULTS 

In this section we explore two different but related issues. 
The first is whether the non-linear behaviour seen in nu- 
merical simulations matches the predictions of analytical 
approaches. The second is whether it is possible to sepa- 
rate the effects of non-linear evolution from primordial non- 
Gaussianity to a sufficient extent for this method to be useful 
in practice. 



4.2 Non-linear evolution and primordial 
non-Gaussianity 

In Fig. |3l we focus on the differences between Avk with 
/nl = 1000 and that with /nl = at z = 0. The per- 
turbative predictions are also plotted for comparison. The 
deviation from the perturbative predictions becomes signif- 
icant as the smoothing scale R is smaller (cro increases) due 
to the primordial non-Gaussian effect coupled with non- 
linear gravity. The increase of deviations at larger cro are 
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V V 

Figure 2. Four Minkowski Functionals Vj. {left) and their difference ratios Aiij. (right) for the simulated mass density fields at z = 3.96, 
with /nl = (filled circles), 1000 (open triangles) and —1000 (crosses). The simulated fields are smoothed with a Gaussian window 
function at the scale R = 10/i~^Mpc. The error-bars denote the sample variance estimated from 1000 Gaussian realizations with same 
z,R and box-size as the simulations. For comparison, the theoretical expectations from perturbation theory (equations [2] and [3) are 

nlnttpH with line's 
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Table 1. Root- mean-square differences between simulated MFs and perturbative formulae at different z and /nl with the corresponding 
(TO- The smoothing scale is fixed at R = 10/i~^Mpc. 





JNL 


^0 


Vb 


Vi 


V2 


V3 


Vb 


Vi 


V2 




5.15 





0.080 


0.003 


0.006 


0.009 


0.024 


0.15 


0.27 


0.24 


0.41 


5.15 


1000 


0.080 


0.005 


0.008 


0.013 


0.029 


0.17 


0.27 


0.24 


0.38 


5.15 


-1000 


0.080 


0.003 


0.006 


0.012 


0.024 


0.33 


0.44 


0.40 


0.54 


3.96 





0.099 


0.003 


0.006 


0.009 


0.023 


0.14 


0.24 


0.20 


0.37 


3.96 


1000 


0.099 


0.006 


0.010 


0.014 


0.031 


0.20 


0.34 


0.26 


0.41 


3.96 


-1000 


0.099 


0.003 


0.004 


0.008 


0.020 


0.14 


0.20 


0.23 


0.41 


2.13 





0.16 


0.007 


0.014 


0.021 


0.044 


0.18 


0.34 


0.29 


0.44 


2.13 


1000 


0.16 


0.011 


0.021 


0.029 


0.057 


0.23 


0.40 


0.34 


0.48 


2.13 


-1000 


0.16 


0.004 


0.010 


0.017 


0.034 


0.15 


0.31 


0.27 


0.40 


0.96 





0.24 


0.015 


0.030 


0.043 


0.081 


0.26 


0.48 


0.41 


0.58 


0.96 


1000 


0.24 


0.020 


0.039 


0.055 


0.10 


0.31 


0.54 


0.48 


0.65 


0.96 


-1000 


0.24 


0.010 


0.023 


0.034 


0.065 


0.22 


0.43 


0.36 


0.51 








0.38 


0.035 


0.067 


0.095 


0.17 


0.40 


0.68 


0.63 


0.86 





1000 


0.38 


0.042 


0.078 


0.11 


0.19 


0.44 


0.73 


0.69 


0.93 





-1000 


0.38 


0.028 


0.058 


0.083 


0.15 


0.36 


0.64 


0.58 


0.80 




Figure 3. The difference of Avf; with /nl = 1000 from those obtained with Gaussian initial conditions At)fe(/NL = 0) at z = for 
different smoothing scales R = 20h~^Mpc (ctq = 0.17), 10/i~^Mpc (ctq = 0.38), and 5/i~^Mpc (erg = 0.74). Simulated results averaged 
over three bins are plotted with symbols and the perturbative formulae are also plotted with lines. 



also seen quantitatively in Table [T] The shape of the devi- 
ation is skewed to the positive side of u with a higher peak 
a.t V — ~l/ao (the number density is zero) while the overall 
amplitude of the deviation Avk is roughly the same as that 
from the perturbative predictions. 

It is interesting to estimate the sensitivity of the MFs 
to primordial non-Gaussianity in the non-linear regime, be- 
cause the effect of primordial non-Gaussianity on the MFs 



should become increasingly significant as the sample vari- 
ance decreases, i.e. at smaller smoothing scales. The MFs 
are, however, strongly correlated with each other among 
different bins of the threshold u and it is therefore neces- 
sary to take into account their covariance when estimating 
the significance of the primordial non-Gaussian effect with, 
e.g., chi-squared statistics. If the covariances among differ- 
ent bins were not considered, one would overestimate the 
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value of chi-square as the total number of bins increases. 
When the field follows nearly Gaussian statistics, the covari- 
ance matrix is well approximated with the one numerically 
estimated from a large number of Gau ssian realizations. 
l|Komatsu et al.ll2003l : lHikage et al.ll2006l ). When the field is 
non-linearly evolved, it is an exceptionally time-consuming 
process to generate enough number of realizations to com- 
pute the inverse matrix of the covariance (the number of 
realizations must be larger than the degree-of-freedom at 
least). 

Instead of calculating the covariance matrix directly, 
therefore, we instead estimate the amount of information 
contained in each MF as a function of v. For this purpose, 
we calculate the effective number of bins N^ff for each MF 
and for all MFs combined as follows: 



(19) 



where i and j denote the binning number of different v and 
different kinds of MFs and A'^bin denotes the total number 
of bins. The covariance matrix dj = (AviAvj) is computed 
from 1000 Gaussian realizations with the same cosmological 
parameters and the same box-size as those of the A*'-body 
simulations. As A'bin is increased in a fixed range of v from 
—3.6 to 3.6, the values of A'eff converges to 2, 6, 8, and 12 for 
each MF from fc = to 3 and then 12 for all MFs combined. 
The results indicate that the correlations among difi'erent 
bins of V is very strong for Vb and that higher fc-th MFs 
have more independent information as a function of v. 

Applying the value of A'eff for non-linearly evolved sim- 
ulations, we calculate the chi-square values of the primordial 
non-Gaussian effect on MFs as a function of /nl as 



. N,a ^ (Aw, 
X (/nl) = ^ 



(SIM) 



(/n 



Av, 



(SIM) 



(/nl = 0)) = 



/A (SIM) /J. 



0)2) 



(20) 



The variance (A«P'^^'(/nl = O)'^) is estimated from 10 real- 
izations of A'-body simulations with Gaussian initial condi- 
tions (the cosmological parameters and simulation box-size 
are the same as for the A''-body simulations). The normal- 
ized MFs AvkifNh) at arbitrary /nl is linearly interpolated 
using the simulation results with /nl ~ and 1000. We 
confirm that the linear interpolation works well using simu- 
lations with |/nl| = 100 and 500. 

Table [2] lists the value of /nl at different R when 
the effect of the primordial non-Gaussianity is compara- 
ble to the sample variance, that is ~ 1- The volume 
of the simulation box-size is 0.125(/^"^Gpc)^ which is less 
than half the volume of the SDSS main galaxy sample 
0.3(/i~^Gpc)'^. As the smoothing scale decreases, the primor- 
dial non-Gaussianity becomes significant. At i? = 5/i~^Mpc, 
the primordial non-Gaussianity with /nl ~ 50 is compa- 
rable to the sample variance and then corresponds to the 
current observational constraints from WMAP. Note that 
the detectability of primordial non-Gaussianity from actual 
observations is, however, strongly dependent on the uncer- 
tainty of the cosmological parameters and the galaxy bias- 
ing, which we have not attempted to model in detail. 



Table 2. The values of /nl at = 1 when the effect of the 
primordial non-Gaussianity is comparable to the sample variance 
of mass density field for different smoothing scale R (equation 
[201 '). The values of N^g for each MF Vk are 2, 6, 8, and 12 from 
fc = to 3 and 12 for all MFs combined. The volume of the density 
field is a cube at a length 0.5/i~^Gpc and the redshift is 0. The 
other cosmological parameters are fixed to be fiducial values. The 
effective number of bins N^g (equation 1191 ) is also listed in last 
line. 



_R[h-iMpc] 






/nl at 


= 1 




Vb 


Vi 


V2 


Vs 


All MFs 


30 


770 


480 


520 


370 


350 


20 


420 


300 


310 


210 


210 


10 


190 


180 


140 


150 


110 


5 


90 


80 


90 


60 


50 



5 SUMMARY AND CONCLUSIONS 

We have studied the imprint of primordial non-Gaussianity 
on the topological properties of LSS using the MFs. Charac- 
terizing primordial non-Gaussianity as a quadratic correc- 
tion to the primordial potential fiuctuation with constant 
amplitude /nl , we compare the MFs with difi'erent values of 
/nl from the mildly to the strongly non-linear regime using 
high-resolution A'^-body simulations. Perturbative formulae 
of the MFs based on the multivariate Edgeworth expansion 
well reproduce the MFs of simulated mass density fields in 
the weakly non-linear regime. When the amplitude of the 
density fiuctuation ctq < 0.1 and |/nl| < 1000, the devia- 
tions of the perturbative formulae from simulations are less 
than a few percent of the amplitude of each MF. They are 
also 10 ~ 40 percent in respect to the non-Gaussian contri- 
butions alone. 

As the fiuctuations become more strongly non-linear, 
the simulated MFs begin to deviate significantly from the 
perturbative predictions owing to non-linear gravitational 
evolution. In order to include small-scale information in re- 
alistic cosmological data sets, detailed numerical analysis is 
therefore essential. 

When we include information from smaller scale fiuc- 
tuations, the effects of primordial non-Gaussianity are in- 
deed significant. Using statistics, we find that the pri- 
mordial non-Gaussianity with /nl ~ 50 has a significance 
level corresponding to la, considering the sample variance 
of mass density fields at _R = 5/i~^Mpc with a volume of 
0.125(/i~^Gpc)^. This implies that measuring the MFs in 
a SDSS-like survey could constrain /nl at a level compa- 
rable with current CMB limits. This is an interesting re- 
sult, since other observations, like the cluster abundance, 
that can effectively constrain /nl at high redshifts, become 
useless at 2 = when non-Gaussian features generated by 
non-linear dynami c s completely obli terate primordial ones 
(jGrossi et al.ll2007l : IXaiig et al.^l2007^ . 

The actual detectability of the primordial non- 
Gaussianity is, however, strongly dependent on the degen- 
eracy between the cosmological parameters and the primor- 
dial non-Gaussian effect. Understanding the properties of 
the galaxy biasing is also very important in determining the 
primordial non-Gaussianity accurately. We will consider this 
issue in a forthcoming paper. 
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